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We present new exact expressions for a class of moments for the geometric Brownian 
motion, in terms of determinants, obtained using a recurrence relation and combinatorial 
arguments for the case of a Ito’s Wiener process. We then apply the obtained exact formulas 
to computing averages of the solution of the logistic stochastic differential equation via a 
series expansion, and compare the results to the solution obtained via Monte Garlo. 

I. INTRODUCTION 

The geometric Brownian motion is the stochastic process described by the differential equation 

df = + afdWt, (1) 

where Wt is a Wiener process and a are constants describing the drift and the variance of the 
noise, respectively. The solution can be written as 

/(IW,t) =exp|^/r-t + (Tlw| . (2) 

Geometric Brownian motion is used for modelling many phenomena in a variety of contexts [1]. 
A prominent role is played in financial applications, where the distribution of returns can be 
approximated by a log-normal distribution [HE], at least in specific regimes. 

For the computation of certain properties, it is necessary to compute the integral of f(W, t) 
over a time interval 

F[W,t]= f f{Ws,s)ds. (3) 

Jo 

The evaluation of this functional is also involved in the solution of the geometric Brownian motion 
with logistic corrections. In general, averages of the form 

OO OO 

{G{F[WM = Y,ak{F[W,tf) 

k=0 k=0 


( 4 ) 
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are quite common. The evaluation of averages of powers of the integrated exponential Brownian 
motion, then, is instrumental for the computation of these observables. 

Detailed studies of this functional and of its powers [an are already available in the literature. 

In this paper, we will derive exact formulas for the evaluation of these integrals, under the 
assumption of the Ito formulation for the Wiener process. Similar results have been given in mm- 
Motivated by obtaining exact formulas for Asian options, in [S] Yor obtained an exact formula in 
terms of polynomials for the following moments: 

rt 




dt 


(5) 


Using Girsanov’s theorem [6], one can derive a series of identities, in which the last is Bougerol’s 
formula 


where 






2n 


( 6 ) 


and 


Pnix) = Tin)Y^ 

j=0 


n 


CjZ^, 




(7) 


( 8 ) 


k^jO<k<n 

In this work, we take a different route with the use of combinatorics. We prove a recurrence 
relation for the integrals involved at the A:-th order in terms of integrals at the {k — l)-th order, 
and after resummation, we get an identity in terms of a determinant. 


II. CALCULATION OF MOMENTS 

The central quantity of interest in the present paper is given by the average over the Wiener 
process Ws'- 


rk{n,a,t) = {F[W,tf). 


(9) 


If we expand Eq. , we obtain 


{F[W,t)= / dik 

Jo 


dti { e 




( 10 ) 
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We will use the following formula due to the properties of integrals with Gaussian measure mm. 
and in which we assume that {WtWf) is of the Ito type. This implies 






( 11 ) 


By using this property, we can now prove the following fact: 

Lemma 1. For the average over the Wiener process Wg of Ito type, the following formula holds 
true: 


t _ pik — 1 

uik / p(m+< 7 ^)T -2 


rk{p,a,t) = r(n) / e^"'* / e 
Jo Jo 

Proof. By a direct application of Eq. (0 


L 


{u+ik-2)P)t2 f ^{fj.+{k-l)P)ii 

Jo 


dii ■ ■ ■ dik (12) 


rk{ia,a,t) = (J ■■■ J e 


■■■dik 


k 

rt rt 


eEti[A-^)h+V E, 


(13) 


JO Jo 

Due to symmetry of integrand, we can order the integration variables as U < ti+i, obtaining 

ft pk-l ^2 


rk{iJ.,a,t) = T{k) 


gEi=i vii+'^2 Ei7tj;i,i=i 


0 Jo 


>0 


(14) 


whence: 


rfc(/i, a, t) = r(A:) f • • • f" 

Jo Jo Jo 

= r{k) f • • • f" j:tiik-i)u]di ^... ^4, (15) 

Jo Jo Jo 

After rearranging carefully the terms, we arrive at the final result: 

rfe(ir, a, t) = r(A:) f ■■■['" e(^+(^-2)-Ai2 

Jo Jo Jo 
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(16) 


□ 


Let us now expand further on Eq. ( |16[ ). It is convenient to hrst perform the rescaling ti = tui. 
Then, by dehning \j = t{p + {k — j)cr^), we obtain 

f'^k r'^2 ^ 

rk{X,t) = t^r(A;) / duk / duk-f- / drtie^»=i = t^r(fc)sfc(Ai, • • • , Afc) (17) 
Jo Jo Jo 
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Therefore, the computation of Q reduces to the computation of 

pi r'^2 ^ 

Sfe(Ai, • • • , Afc) = / duk duk-i--- 


(18) 


"l fUk ru2 

duk-i■■■ 

Jo Jo Jo 

a very similar formula had been obtained in [3]; the main difference between the treatment made 


there and our relies on what follows. An important observation will enable us to evaluate these 
integrals exactly by means of combinatorics: it is possible in fact to prove the following result, 
which establishes a recursion relation among the s^- 


Lemma 2. For the quantity Sk{\i, ■ ■ ■ , A^), we have 

• • • , Afc i) — Sfc i(Ai, • • • Xk-2, Xk-i + Xk) 


'SA:(Ai, ■ ■ ■ , Xk) — 


Xk 


(19) 


with So = 1- 


Proof. Let us write 


pi r'^k 

* * * ? ^k) / duf^ I 

Jo Jo 

= [ dttfce^''“'=/(uA:), 

Jo 


r'^2 


duie^i=^ 


( 20 ) 


with f{uk) = duk-i ■ ■ ■ Jq^ duie^i=i Integrating by parts. 


•Sfc(Ai, • • • , Xk) — 




Xk 


-fiuk) 


^ /*! 

- / duk— — /(uk), 
0 Jo '^k 


( 21 ) 


where 


f{Uk) = 


C'^k-l 


duie'^i^^ AiMi+Aj,_i 


'^k 


— gA(j_iUj. 


JO 

ruk-i 


r“2 


duie^^= 


Jo Jo 

Finally, the identity /(O) = 0 gives the desired result. 


( 22 ) 

□ 


Eq. (19) suggests the evaluation of the averages by means of combinatorial considerations. 


Indeed, the evaluation of the integral can proceed graphically, for any fixed order of the moment 
k, as in Fig. Starting from the top and using the properties of the recurrence relation, at each 
order k, we organize the recurrence on a binary tree. In Fig. each left branch will pull out a 
factor e^/A, where A is the A obtained from the previous order. One has to consider the fact that 
however, in each right branch one pulls out a factor — 1/A, and sums the two factors of A’s. To 
obtain the final formula, once the empty set has been reached, one multiplies the final term by all 
the factors in the branch. 







5 


o o 


—O O Q 1 

^^2 A2 + A3 ^2 + A3 

o 00 

>» •>» 


gAi gAi+As 1 1 gAi+A2+A3 I 

Ai A] Ai 4- A 2 Ai + A 2 Ai A] Aj + A 2 + A 3 Ai + A 2 + A 3 

FIG. 1: Application of the graphical method for evaluating k-th moments of the exponential integrated 
Gaussian process to the case k = 3. 


We now give exact formulas for all the terms obtained from the recurrence. To lighten the 
notation, we define Hkj = Xj + ■ ■ ■ + Xk and Sk = Sfc(Ai, ..., A^). 

By iterating the second term of the recurrence, we obtain 

Sk = - Sk-l - Sk-2 H-'Sfc_2(Al, . . . , Xk-3, fJ-k,k-2) 

f^k,kf^k^k—l l^k^kl^k^k—1 

gMfc,/c 1 gMfcjfc—2 

= - Sk-l - Sk-2 H- Sk-3 -Sfc_3(Al, . . . , Afc_4, Hk,k-3) 

t^k^k f^k,kf^k^k—l f^k^kf^k,k—ll^k,k—2 f^k,kl^k,k—lf^k,k—2 


—^,- 1 +(-1)^—^—. 

^^k,j k'kj+l k‘k,k Mfc,l l^k,k 


Therefore, 


oAjH-hAj. 


■SfcCAi,..., Xk) = ^ ^(—1)^ ^—r-So—i(Ai,..., Ao—i) H-r-^^-. (23) 

^ ic ’ + + ^ nti(A. + --- + A,) 


Lemma 3. Let fk = Yl^=o^j '^hen 


m=0 0=io<ii< — <im<k+l=im+i 


g^m+l Ig^m 1 , , , 1 

—1 ^0 


Proof. We proceed the proof by induction on k. We have fk = Cq, which holds for k = 0. We 
assume that the claim holds for k and let us prove it for k + 1. By the recurrence relation, we have 


/fc+i = e^+^ + ^e^^+Vj_i. 

i=i 
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By induction hypothesis, we have 


k-\-l 


j-1 




e| 'ej'" '■■■£*; ‘ 


= e, 


fc+i 


i=i 
fc+i /i-i 

+E E 


m=0 0=io<ii<---<im<i=*m+l 


E 


gim+2-lgi-lgi^-l _ _ _ gil-1 


j=l \7n=0 0=io<ii<---<im<j=irrL + i<im+2=k-\-2 

fc +1 / fc +1 


= a‘+‘ + E e 




—1 ^0 


™=1 \i=m 0=io<h<---<im-l<i=*m<*m + l=fc+2 
fc+1 


= e, 


fc+1 


+ E E 


g^m+1 Igim 1 , , , 1 

—1 ^0 


?T2.— 1 0 — 2Q ^^^TTl ^^771 + 1—fcH”2 

fc +1 

= E E 

m=0 0=io<ii<"-<im<im+l=k+2 

which completes the induction step. 


*m + l-l im-l . . . ph-1 
im im — 1 ^0 ’ 


If e^ = (-i)k-J^ 


.E+'"+^fc 1 i. (-i)*^ 


and Cn = 


n+(^iH—1"^*:) ^ n^=i(EH—hAfc) 

fc 


then Eq. (23) can be written as 


Sfc(Ai,..., Afe) = ^ ejSj-i(Ai,..., Aj-i) + eg. 
i=i 

Applying Lemma we obtain 
Sfc(Ai, • • •, Afc) — 


□ 


S S 

LIT.—0 0—2q ^2 ^ ^ fc+1—2^_|_ 


g^TTi-l-l Ig^TTi 1 ^ ^ ^ g2l 1 

—1 ^0 


(- 1 ) 


21-1 


<fc+i=i^, 1 n}Li^(Aj H h Aq_i) ]Tj=i^ ^(Af ^-h Ai^_^j_i) 


n 




hAi^ 1 


rh+i —1 / 


= E E 

m=0 0=io<n<-"<*n 
k 

= E E 

m =0 0 =io<n<-<im<fc+l=im + l n}L/(Aj H-h Aii_i) n^i (Aj -h Ai^_^j_l) 

which leads to the following result. 


-hAfc 


Theorem 4. For all k > 0, Sfc(Ai,..., A^) is given by 


E E 


-l-Aj, 


m=0 0=io<ii<-<im<k+l=i„^+l rijLl (Aj H-h Aii_i) n^l Y\j=it (Aj H-h Aj^^^_i) 

The above theorem gives 


^fc(^Ai,..., tAfc) 

k 

E E 

m=0 0=io<ii<-"<im.<k+l=im+i 


^_-['jk-m^t(Xi^A -hAfc) 


+ ■ ■ ■+A..-,) nr.i ni;y‘(A, +... + ' 
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Hence, the coefficient of t” in Sk{tXi ,..., tXk) is given by 
{n + k — 1)! ^ 


E 


m=0 0=io<ii<-<im<k+l=im+l H-^ ^il-l) WT=l H-^ ^i«+l-l) 

A cleaner formula can be obtained in terms of determinants. If we define Sk = Sfc(Ai,..., Xk) for 


all /c > 0, where we define sq = S-i = 1. If = (—1)^ ^ 


AjH- \-\f. 


e J 


- 


and Cq = 


- ’ 


then Eq. (23) can be written as 


j=0 

Using Theorem 4.20 in [8] (see also M), we obtain the following result: 
Theorem 5. For all k > 0, = Sk{Xi ,..., Xk) is given by 


( 


det 


^k-l ^k-2 
^ ^k-l ^k-2 


\ 


0 -1 

0 0 

0 0 


pfc-2 

-A:-2 


^k—1 ^k—1 ^k—1 


^k—2 ^k—2 ^k—2 


where e^ = (—1)^ ^ 


AH-l-Aj, 


and Cn = 


-1 el 
0 -1 
(- 1 )" 


(24) 


So / 


nL,(Ai+-+AU f3fe^^(A,+...+Afe)' 


This expression can be elucidated by means of an example. Consider the case n = 3 and a = 0. 
In this case, we have an exact formula for Xj = fit: 


n = 3 : Det 


/ 

1 

•>-. Irt> 

^ 1 : 

- 1 / 2 ^ 

1/6 

“1/6 

\ 


-1 

eE: 

fl t 


1/2 7^ 



0 

-1 


_ j_ 



/i t 

(it 


V 

0 

0 

-1 

1 

/ 


= 1/6 


(eMt _ 

fFD 


which provides the exact value obtained from the deterministic logistic equation. 
The previous results are general enough to hold also for averages of the form 

hn{fi,a,t) = )*(e'"^‘T[IT, t]"). 


(25) 


It is easy to see that Eq. (24) applies. We just substitute = Afc + ta^ and multiply by a factor: 


rn{fi,(T,t) = e^^ "2 *r(n)s„(Ai, • • • , An) 

= e^*rr(n)sn(Ai, • • • , An). 


( 26 ) 
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FIG. 2: Plot of the functions s„(/i, cr, t) and s„(^, a, t) provided in Appendix [A] as a function of t ioi fj, = 3, 
(T = 1 for n = 1, 2, 3,4. 


The above expression will be used in the following section as an application to averages in the 
logistic stochastic differential equation. Appendix [^contains the analytical values of the functions 
■Sfc(Ai, • • • , Afc) and Sfc(Ai, • • • , A^) for the cases n < 4. It is immediate that, when t 0, we have 

+ cr^), • • • , t(/i + (n - 1)(T^)) (27) 


This facts implies that Eq. (26) is analytic in t 


0. Fig. [^is a plot of Sn{t) for /r = 1, u = 0.1. 


III. AVERAGES OF LOGISTIC SDE IN PERTURBATION THEORY 


As an application of the formula (26), we focus on the solution of the logistic stochastic differ¬ 
ential equation motion, 


dx = X 


given by (we follow mm) 


— — ]dt + adW 


-1 


:{t) = XQ ^ y 

„-n \ ^ Jo J 


We evaluate the average of the solution x{t), 


(28) 


(29) 


{x{t)) = ( xo ^ J 


= Xo 


E(-i) 


n=0 


n=0 


\ X / 


= Xo e 


fit 


oo 

^(_1)- ^t^yr{n)sn{Xu--- ,An). 


n=0 


(30) 
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FIG. 3: Plot of the mean of solution of the logistic stochastic differential equations for /r = 1, tr = 0.1, 
X = 100, and with dt = 10“^, solved numerically using a stochastic Euler method and averaged over 
1000 simulations (solid red), versus the analytical solution obtained at fc-th order, considering ^ ^ the 

perturbative parameter. We can observe that at higher order we obtain a solution closer to the one simulated. 


and observe that it involves the moments of (26). In the limit x —?• oo, Eq. (30) reduces to the 
average of the geometric Brownian motion. We consider now truncations of the mean of x{t) 
at A:-th order, {x{t))k = xq ‘ and compare the truncated 

solution to the one obtained numerically. In Fig. we plot {x{t))k for A: = 0,1,..., 4 obtained for 
H = I, (7 = 0.1, X = 100 by means of a stochastic Euler method with dt = 10“^, and the averages 
then obtained using Monte Carlo over 1000 samples. We observe that the higher the order of the 
approximation, the closer we are to the solution obtained by Monte Carlo. 


IV. CONCLUSIONS 

In this paper, we have presented new exact formulas for the moments of the integrated expo¬ 
nential Brownian motion in terms of sums and determinants, and based on recent results obtained 
in [9j. We described a simple graphical method to evaluate them, based on a recurrence relation. 

Exact formulas were proved in [5] in terms of polynomials. In this paper however, we have 
taken an alternative route based on combinatorics. After realizing that the mean can be evaluated 
exactly using the properties of Gaussian integrals, and after observing that these moments feature 
a recurrence relation, we have shown that exact expressions can be obtained via a combinatorial ar- 
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gument. These exact expressions were then observe to be equivalent to evaluating the determinant 
of a specific linear operator which depends on the order of the moment. 

To complete the presentation, we have applied the formulas to the exact solution of the logis¬ 
tic stochastic differential equation. There, the evaluation of the ensemble expectation values of 
certain observables can be carried out with our method, via Taylor expansion. In particular, the 
comparison of the mean solution obtained by means of Monte Carlo simulations ^ with our method 
shows that our formula permits to approximate to a higher precision some ensemble averages of 
properties of the solution of the stochastic differential equation. 
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Appendix A: Exact formulas for n < 4 
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